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ABSTRACT 



We present the result of an indirect search for high energy neutrinos from Weakly 
Interacting Massive Particle (WIMP) annihilation in the Sun using upward-going muon 
(upmu) events at Super-Kamiokande. Data sets from SKI-SKIII (3109.6 days) were used 
for the analysis. We looked for an excess of neutrino signal from the Sun as compared 
with the expected atmospheric neutrino background in three upmu categories: stopping, 
non-showering, and showering. No significant excess was observed. The 90% C.L. upper 
limits of upmu flux induced by WIMPs of 100 GeV c~^ were 6.4x10"^^ cm~^ and 
4.0x10"^^ cm~^ for the soft and hard annihilation channels, respectively. These 
limits correspond to upper limits of 4.5xl0~^^ cm~^ and 2.7x10"^'' cm~^ for spin- 
dependent WIMP-nucleon scattering cross sections in the soft and hard annihilation 
channels, respectively. 

Subject headings: dark matter, — neutrinos 



Introduction 



From recent observations, Weakly Inte racting Massive Par ticles (WIMPs) are considered a 



favorite candidate for cold dark matter (jBennett et al.l |2003| ). From the viewpoint of mini- 



mum supersymmetric extensions of the Standard Model, the most well-motivate d candidate for 



WIMPs in the univer se is the lightest supersymmetric neutral particle (LSP) (jJungman et al 



1996 



BergstromI l200dl . One c andidate for the LSP which is theoretically well-studied is the neu- 



tralino (x) (jPrees et al.lll993l ). The neutralino is a linear combination of supersymmetric particles 
that mix after electroweak-symmetry breaking. Taking into account constraints from accelerator 
experiments and astronomical observations, the neutralino mass is expected to be in the range from 



— ^" • 1 

(Amsler et al. 


2008: 


Gilmore 


2007) 



One method of searching for a WIMP dark matter signal is an indirect search where decay or 
annihilation products from WIMPs are observed as originating from the center of a gravitational 
potential well such as a celestial body. WIMPs have a small but finite probability of elastic scatter- 
ing with a nucleus. As matter is concentrated in deep gravitational potential wells, this scattering 
is more likely to occur near the centers of such wells; if their final velocities after the scattering are 
less than the escape velocity, they are trapped gravitationally and eventually accumulated in the 
center of the gravitational potential. There are two types of scattering: spin-dependent (SD) and 
spin-independent (SI) scattering. In the former, WIMPs couple to the spin of the target nucleus, 
and in the latter WIMPs coherently scatter on the nucleus with probability scaling as the square 



Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 



- 4 - 



of the mass of the target nucleus. WIMPs, being Majorana particles, can annihilate in pairs, and 
produce primarily r leptons, 6, c and t quarks, gauge bosons, and Higgs bosons depending upon 
their masses and compositions. At the final step, the decay of the WIMP annihilation products 
may produce many kinds of particles such as neutrinos, positrons, antiprotons, antideuterons, etc. 

In the Sun, WIMPs are most likely to scatter on the dominant component, hydrogen, either 
via SD or SI coupling. Direct detection searches on Earth turn out t o have a sensitivity a dvantage, 
generally, for SI interactions, while indirect searches apply to both (jKamionkowskil Il995l ) . WIMP 
annihilation products in the form of muon-neutrinos are an excellent instrument for indirect searches 
since they can pass through the matter of the Sun, and interact in the Earth; the resulting detected 
muon will be correlated with its parent neutrino direction, providing an unmistakable signal. In 
the past, several neutrino telescopes such as Super-Kamiok ande, AMANDA, and IceCube have 



reported the result s of indirect searches for WIMPs in the Sun (jDesai et al.ll2004l : iBraum et al. 



200S 



Abbasi et al.ll2009l ). Here we investigate upward-going muons (upmus) which are generated from 
high energy neutrinos come from the WIMP annihilations using the Super-Kamiokande (Super-K, 
SK). We search in the direction of the Sun, and an excess of neutrino flux above the atmospheric 
neutrino background is sought in the upmu events. Although recent new detectors such as IceCube 
have larger acceptances for higher mass WIMPs, SK is better equipped for the search of lower mass 
WIMPs (<100 GeV) due to its lower energy threshold for neutrino signals. 

A previous WI MP search in SK was done in 2004 using a dataset of 1679.6 days of through- 
going upmu events (jDesai et al.ll2004l ). In this paper, we have nearly doubled the size of the dataset 
(3109.6 days). In addition, we have divided the upmu events into three categories, taking into ac- 
count their neutrino energy dependence of detector response in order to improve the limit on the 
WIMP-induced upmu flu x and SD cross secti on, especially for lower mass WIMPs. This analysis 
also now uses WIMPsim (jeiennow et ahlboosl ) and the DARKSUSY (jOondolo et al.ll2004l ) simula- 
tor in the calculation in order to simulate the WIMP-induced neutrino fluxes and propagation. 



2. Indirect WIMP search in Super-Kamiokande 

The SK detector is a 50-kton water Cherenkov detector located at the Kamioka mine in Japan. 
Since SK started operation in 1996, there have been four experimental phases: SK-I (1996 April- 
2001 July), SK-II (2002 October-2005 October), SK-III (2006 July-2008 August), and SK-IV (2008 
September-). We use the datasets from SK-I to SK-III (3109.6 days) for this analysis. 

Upmu events, which are categorized as the most energetic events in SK, are used for the 
analysis. Muon neutrinos interact in the rock surrounding the detector to produce muons via 
charged current interactions. If the energy of the incoming neutrino is high enough, the muon 
produced by the neutrino will enter the detector and emit Cherenkov light. To separate neutrino- 
induced muons from cosmic ray muons, we select only upmus, since the background from downward- 
going cosmic ray muons overwhelms any neutrino-induced muons from above. Both the neutrino 
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cross section and the range of the produced muon are proportional to the neutrino energy. Thus, the 
effective target volume is very much increased if we use the upmu events for high energy neutrino 
searches. Upmu events are divided into three categories: stopping, non-showering through-going, 
and showering through- going. If upmus stop in the detector, they are categorized as "stopping." 
Stopping upmus are the lowest energy portion of the upmu sample. Upmus which fully traverse the 
detector are classified as "through-going." The category of through-going is then further subdivided 
into "showering" (if the event has accompanying radiation due to radiative energy loss) and "non- 
showering" (other wise). The detai ls of the showering versus non-showering classification scheme 
are described in (jPesai et al.ll2008l ). Showering through-going events make up the highest energy 
events in the upmu sample. The detector angular resolution is 1° for stopping and non-showering 
through-going, and 1°.4 for showering through-going muons, as estimated by Monte Carlo (MC) 
simulation. This analysis is based on a total of 4351 upmus: 919 stopping, 2901 non-showering 
through-going, and 531 showering through- going. 

The main background source for this analysis is upmu events produced by atinosphe ric neu- 
trinos. These were simulated by the NEUT neutrino simulation package iHayatol EqO^) and a 
GEANT-based detect or simula t or. W e use the neutrino flux calculated by ( iHonda et al.l (j2007l ) 
; El, <10 TeV) and ( IVolkoval (jl980l ) ; Ei, >10 TeV) for this simulation. The main uncertainty 
in the absolute number of upmus comes from the uncertainty in the cosmic ray flux (~ 20%). In 
this analysis, this uncertainty is irrelevant since we normalize the absolute number of simulated 
atmospheric neutrino events by the number of upmu data events in each category. The effect of neu- 
trino oscillation is also taken into account in the simulation. We consider — t^t oscil lations with 



Am^ = 0.0025 eV^ and sin ^ = 1.0 independently obtained from the previous analyses (jAshie et al 
200,5l . l2004l ). 



3. Results of the WIMP search 

3.1. Upmu events from the Sun 

The event distribution (SKI-III:3109.6 days) as a function of cos^sun for each category of 
upmu is shown in Figured) Here, 9 is the angle between the reconstructed upmu direction and the 
direction of the Sun. This 9 is called the "cone half-angle." In Figure [H cos Sgun = 1 is defined as 
the direction of the Sun. To compare data with the MC events, each MC event was assigned a time 
sampled from the observed upmu arrival time in order to match the livetime distribution and the 
detector acceptance of the observed events. The number of MC events is normalized to the total 
number of data events for each category. Here the normalization was done for the entire cos^sun 
region for simplicity. Without signal regions the difference of normalization is less than 1%. 

We search for an excess from the direction of the Sun in each of the three upmu categories: 
stopping (A), non-showering (B), and showering (C). The fraction of each upmu category as a 
function of parent neutrino energy {Ejj) is shown in Figure [2j For typical energies of 10, 100, and 
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1000 GeV, they are: A:B:C = 0.73:0.20:0.07, 0.10:0.62:0.28, and 0.02:0.50:0.48, respectively. 
These fractions were obtained from the atmospheric neutrino MC simulation. The fraction depends 
only on parent neutrino energy; its source is irrelevant. We can roughly deduce the WIMP mass 
from the event category where an excess is observed. Since the spectrum of the atmospheric 
neutrino background is proportional to E'"^'^, while the signal spectrum from WIMPs gives a 
bump structure, this energy-dependent approach is good for reducing the atmospheric neutrino 
background. We separately analyze the three upmu samples, then combine the results by taking 
into account each fraction assuming a WIMP mass. 

Since the directional correlation between neutrinos and their produced leptons depends on 
neutrino energy, it also depends on the WIMP mass and its decay mode. Therefore we look for a 
possible correlation of upmu direction to the Sun by searching for an excess of events in various cone 
half-angles up to 30°. Using DARKSUSY, we checked the neutrino-muon angular correlation for 
the neutrino spectrum generated by WIMP annihilation in the Sun for various SUSY parameters 
and found that the angle between the two is smaller than 30° for most cases. The number of events 
observed in each cone half-angle toward the Sun is shown in Figure El Zero degrees corresponds to 
the direction of the Sun. Crosses indicate the data from SKI-III (3109.6 days) and the histogram 
indicates atmospheric neutrino MC with same normalization as in Figured) 



3.2. Upmu flux limit from the Sun 

We calculate the upper limit of upmu flux from the direction of the Sun with respect to 
the WIMP mass. Here we estimate the cone half-angles which contain more than 90% of the 
total nui nber of upmus expe cted from WIMP annihilations in the Sun. We use the WIMPsim 



package (iBlennow et al.ll2008l ). which simulates WIMP annihilations, propagation of neutrinos from 
those annihilations, and the propagation of muons after neutrino interactions in the rock. In order 
to derive the expected neutrino spectrum from the WIMP annihilation, we consider two cases: 
annihilation into bb (soft channel), and annihilation into W~^W~ (hard channel). These give the 
softest and hardest neutrino spectra, respectively. The cone half-angles obtained from these spectra 
are shown in Figured! The cone half-angle for the soft channel is generally wider at the same WIMP 
mass since many neutrinos from annihilation via the soft channel are expected to have lower energy. 

The other smearing effects for angular correlation are muon propagation in the surrounding 
rocks, distribution of WIMPs in the Sun, and angular resolution of upmu reconstruction. The 
smearing due to Coulomb multiple scattering in the surrounding rocks is 2° - 3° for muon s with 



energy above 1 GeV. The angular size of the Sun at the Earth is 0°.5, which is negligible (jEdsio 



19971 ). These two smearing effects were taken into account by WIMPsim. The detector angular 
resolution of 1° (stopping and non-showering) and 1°.4 (showering) are smaller than the smallest 
cone half angle and thus were not taken into account. 



We show the cone half-angle for several WIMP masses in Tables [D and [2j Table [2] begins at 80.3 



-7- 



GeV since WIMPs with mass lower than 80.3 GeV cannot annihilate into W^W . We followed 



the DarkSUSY for the treatment of masses of standard model particles (jBlennow et al.l 120081 ) 



For each cone half-angle, which corresponds to a certain WIMP mass in each annihilation 
channel, we calculated the upper Poissonian limit of the upmu flux: 



C.L.l 



e X A X T 



(1) 



where A is the detector cross-sectional area in the direction of the expected signal averaged over Sun 
directions. For each position of the Sun, this area is calculated using a minimal track length of 7m in 
the inner detector as a criterion to select viable upmu trajectories. It varies from 968 m^ to 1290 m 
with the position of the Sun. T is the experimental live time and e is the detection efficiency. For 
upmu samples, the inefficiency of detection is negligible (jPesai et al.ll2008l ). The average exposure 



^ X T is 1.67 X 10^^ m'^s. Ni^q is the upper Poissonian limit on the number of upmu events at 90% 
C.L., given the number of observed events of each category (n^^g), expected background of each 
category (?^BQ) and fraction of each category (F*) where index i corresponds three types of upmu 
category. This fraction is estimated from the neutrino spectrum of the given WIMP mass for each 
annihilation channel and the fraction of the upmu categories as a function of the parent neutrino 
energy as shown in Figure [2j The neutrino spectra were simulated by WIMPsim. Figure [5] shows 
the simulated neutrino spectra for the soft and hard annihilation channels, assuming a neutralino 
mass of 100 GeV. 

iVfln was calculated combining all categories of upmu for a certain mass based on Bayes' theo- 
rem (lAmsler et al.ll2008l ). The Poisson-based likelihood function for the expected signal is as follows 



where Vg is the expected real signal. Using the likelihood function in Equation ( [2]), the upper 
confidence level is written as 

fc"LKbJ^s)dz^s 
L,=oL(nobsl'^s)di/s 

Applying 90% as the C.L. in Equation ([3|), the upper Poissonian limit (A^'go) is obtained. 

The summaries of the numbers of upmu events and the estimated backgrounds for the cone 
half-angles which correspond to each WIMP mass, A'^go, and the calculated flux limits are shown 
in Tables Hand [21 
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The 90% confidence level flux limits as a function of WIMP mass are show n in Figures El and [71 
Results from other experiments, AMANDA (|Braum et al.ll2009l ) and IceCube (lAbbasi et al.ll2009l ). 
are also shown in these figures. The previous Super-K result (|Desai et al.ll2004l ) is also shown in 
Figure [U but please note that it was calculated assuming an annihilation branching ratio. That 
result assumed that 80% of the annihilation products were bb, 10% were cc, and 10% are rf, as 
a typical branching ratio. We cannot regard this as a pure soft annihilation, but it is shown for 
reference. 



3.3. Limit on WIMP-proton SD cross section 



From the upmu flux limits shown in Figures [6] and O we also calculate the upper limit on the 
SD scattering cross section of WIMPs. If we assume that the capture rate of WIMPs and the WIMP 
annihilation rate are in equilibrium in the Sun, the neutrino flux from annihilation depends only 
on the capture rate. Since the age of the Sun is much longer than the timescale of the annihilation, 
it is a safe assumption that they are in equilibrium. We can relate the upmu flux, to the SD 
cross section, o"^^, on the assumption that only one annihilation c hannel is dominant an d the SI 



cross section is zero. The SD cross section, cj^^, can be written as (jWikstrom et al.ll2009l ) 



(4) 



where Kj^(mi^) is the conversion factor between the SD cross section and the muon flux. This value 
depe nds on the WIMP mass, m^. We use the factor introduced in Figure 3 of (jWikstrom et al 
2OO9I ). which was calculated using DARKSUSY. 



The calculated upper limit is shown in Fig ure El Also shown are results from direct det ection 
experiment s: DAMA /LIBRA allow ed region (ISavage et al.l bood ). KIMS (iLee et al.1 boOTI). and 
PICASSO (lArchambault et ahlbood ): and from neutrino telescopes: AMANDA (jBraum et ahlbood ) 
and IceCube (j Abbasi et al.l l2009l ) . All of the limits fr om the neutrino teles copes (including our 
limit) are calculated by the same method described in (jWikstrom et al.ll2009l ). The previous limit 
from Super-K is also shown for reference. Note that the gaps at WIMP masses of ~ 80 GeV 
and 174 GeV in this previous result are due to the thresholds of the W^W~ and rr annihilation 
channels. These do not appear in the new result because we assume only one annihilation channel. 
Our limit is improved at low WIMP mass (below 100 GeV), which is stronger than limits from 
other indirect-detection experiments. Our limit is also stronger than direct-detection experiments. 



4. Conclusion 



We have performed an updated indirect search for WIMPs using the Super-Kamiokande data 
of SKI-III (3109.6 days). A WIMP-induced neutrino signal was searched for in the direction of the 
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Sun. 

We compared upmu events with the estimated background of atmospheric neutrinos. No 
significant excess was observed. 

To calculate the upmu flux limit, we used a new method involving WIMPsim and the DARK- 
SUSY simulator. Cone half-angles which contained more than 90% upmu signal for each WIMP 
mass were extracted from simulation. Upper upmu flux limits as a function of WIMP mass for 
the soft and the hard annihilation channels were obtained. Our limit extends into the region of 
expected upmu flux estimated by some SUSY models. 

An upper limit on the SD scattering cros s section was also obta ined. We used a WIMP mass- 



dependent conversion factor as described in (jWikstrom et al.ll2009l ). At low WIMP mass (below 



100 GeV) our limit on the WIMP-proton SD cross section is better than other experiments for both 
direct-detection and indirect-detection experiments. 
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Fig. 1. — Distribution for cos 0sun of upmu events relative to the Sun. Here, cos^sun = 1 is the 
direction of the Sun. Crosses indicate the data from SKI-III (3109.6 days) and the histogram 
indicates atmospheric neutrino MC with — Vr oscillations normalized by the total number of 
data events in each category. 
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Fig. 2. — Fraction of each upmu category as a function of the parent neutrino energy (Ev\ Three 
upmu categories are shown as fohows : stopping (black sohd hne) , non-showering (red dashed hne) , 
and showering (blue dotted line). The total of the fractions for the three categories is normalized 
to 1. A Monte Carlo simulation for atmospheric neutrinos is used to calculate this. 
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responds to the direction of the Sun. Crosses indicate the data from SKI-III (3109.6 days) and 
the histogram indicates atmospheric neutrino MC with z/^ — i^r oscillations normalized to the total 
number of data events in each category. 
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Fig. 4. — Cone half-angle which contains more than 90% of the total number of upmus expected 
from WIMP annihilations in the Sun when we assume the soft annihilation channel (black solid 
line) and hard annihilation channel (red dashed line). This result was obtained from WIMPsim. 
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Fig. 5. — Simulated neutrino flux (t'^ + z^) at the surface of the detector from the annihilation 
of neutralinos of 100 GcV in the soft annihilation channel alone (left panel) and hard annihilation 
channel alone (right panel). 
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Fig. 6. — Upm u flux limit vs. W IMP mass for the soft annihilatio n channel. T h e flux limits 
of AMANDA ( braum et aP ^2Qod ). ; green dotted line ), IceCube ( Ubbasi et all (12003),; blue 
dot-dash l ine), and th i s anal ysis (red solid line) are shown. Also shown is the previous limit from 
Super-K ( iDesai et al.l (|2004l ). ; black dashed line). Note that SK's previous limit does not assume 
the pure soft channel, but assumes a typical annihilation channel in which the branching ratio for 
the soft annihilation channel is relatively high (80%). The shaded region is the expected flux region 
from DARKSUSY. 
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Fig. 7. — Upm u flux limit v s . WI MP mass for the hard annihilati on channel. The flu x limits 
of AMANDA ( iBraum et all kood ). ; green dotted line), IceCube ( kbbasi et all kood ). ; blue 
dot-dash line), and this analysis (red solid line) are shown. The shaded region is the expected flux 
region from DARKSUSY. 
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Fig. 8. — Limit on the WIMP-proton spin-dependent cross section as a fu nction of WIMP ma ss 
Limits from direct detection experiments: DAMA/LIBRA allowed region ( I Savage et all hood). ■ 



dark r ed and light red filled, for with and with out ion channeling, r espec tively), KIMS ( iLee et al 
(120071 ) ■ ; hght blue crosses), and PICASSO ( Urchambault et alj pOod ). ; grey dotted line) are 
s hown. Also w 'e sho w here the results of indirect detection (neutrino telescopes): AMANDA 
( iBraum et al.l (|2009l ) , ; black line with triangles) , IceCube ( IWikstrom et al.l (|2009l ) , ; blue hue 



with squares), and this analysis (red line with stars). Also the previous limit from Super-K (green 
dashed line) is shown. This plot was made from SUSY Dark Matter /Interactive Direct Detection 
Plotter, http:/ / dmtools.berkeley.edu/limitplots, . 
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Table 1: Event summary to calculate the flux limit {(p) assuming the soft annihilation channel is 
dominant. The numbers of events (nobs) and estimated background (nBc) for each cone half-angle 
(9) which corresponds to a certain WIMP mass (m-x) shown for all categories of upmu. This 
cone half-angle is defined by the criterion that more than 90% of the WIMP-induced upmus are 
contained for the assumed WIMP masses. Also we show the estimated fraction {F) of each upmu 
category for the given WIMP mass. The calculated result of A'go and the flux limit for each WIMP 
mass are shown in the last two columns. 
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Table 2: Event summary and upmu flux limit assuming the hard annihilation channel is dominant. 
The meaning of each column is the same as in Table. [TJ 
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